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The hard contributions to the heavy quarkonium-nucleon cross sections 

are calculated based on the QCD factorization theorem and the nonrelativis- 
tic quarkonium model. We evaluate the nonperturbative part of these cross 
sections which dominates at y / ijv r /v ~ 20 GeV at the Cern Super Proton Syn- 
chrotron (SPS) and becomes a correction at ^Jsnn ~ 6 TeV at the CERN 
Large Hadron Collider (LHC). J/^ production at the CERN SPS is well de- 
scribed by hard QCD, when the larger absorption cross sections of the x 
states predicted by QCD are taken into account. We predict an ^-dependent 
polarization of the x states. The expansion of small wave packets is discussed. 
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The aim of the present paper is to investigate color coherent phenomena in the prop- 
agation of hadronic states with hidden charm in pA and AB collisions. We restrict the 
consideration to the central rapidity region of nuclear beams at yfs w 20 AGeV (CERN 
SPS) where formation time effects in J /^/-production are corrections. Qualitatively new ef- 
fects for v/i « 200 AGeV (RHIC) and » 6 ATeV (LHC) energies will be discussed also. 
We assume that hadronic states with hidden flavor are produced predominantly in hard 
parton collisions and calculate the absorption factor S of the total yield of heavy quarkonia 
which is not distorted by the initial state interaction, cf. jlj. 

Under these assumptions S can be directly measured by the ratio of the yields of heavy 
quarkonium to the dimuon yield 0. We want to emphasize that the agreement of our 
calculations of the J/\I/ production in pA and in AB collisions with the data strongly supports 
the idea that hidden charm states are perturbatively produced in NN collisions. 

The QCD factorization theorem is used to evaluate the PQCD cross sections of heavy 
quarkonium interactions with ordinary hadrons. However, the charmonium states (here 
denoted X) are not sufficiently small to ignore nonperturbative QCD physics. Thus, we 
evaluate the nonperturbative QCD contribution to the cross sections of charmonium-nucleon 
interaction by using an interpolation between known cross sections. The J/ty-N cross section 
evaluated in this paper is in reasonable agreement with SLAC data ||. 

PQCD predictions for the cross section of T-N scattering should be more reliable since 
(differenctly from the J/^/'s) the spatial size of T's is significantly smaller than the typical 
size for soft hadron interactions. 

These X-N cross sections are used to evaluate nuclear effects for the production of X 
states. We find that the suppression of the J/\l> in pA and AB collisions is reasonably well 
described within the present approach. 

We predict a new QCD effect - the polarization of produced Xc states: cc-states with 
nonzero orbital momentum will be polarized due to the longitudinally directed nuclear (color 
filter) target. Their polarization is evaluated quantitatively. Random interactions with 
comovers will depolarize the x states. Thus, measuring the x polarization may turn out to 

2 



be an effective method to investigate their interactions with comovers. 

The production of charmonium states in AB collisions has been proposed (via J/^ 
suppression) as signal for Debye screening in deconfined matter [|J. However, secondary 
collisions of J/^'s with hadrons have been suggested early on as a competing origin of J/^ 
suppression [§,0. This topic has been considered also in basing on QCD motivated 

models (but not on the QCD factorization theorem as in our paper). In j7j (but not in |§) it is 
also concluded that PQCD predicts a too small J/^f-N cross section to explain the observed 
absorption of J/\l/ at CERN energies. In contrast to the "preresonance" hypothesis suggested 
in [0 and to the model in || we explain the J/\l/ suppression observed in AB collisions as a 
consequence of the predicted large absorption of x states. We found that at RHIC energies 
perturbative and nonperturbative QCD contributions to the J/\I/-iV cross section become 
comparable in variance to 0. At LHC energies the hard contribution will dominate (we 
ignore in this paper the unitarity corrections which will be important at LHC, cf. [[IT)| , |TT|j ) . 
The T-7V" interaction cross section is calculated for the first time in this paper. 

The first evaluations of <J to t{J /^f-N) have been obtained by applying the Vector Dom- 
inance Model to J/ 1 ]/ photoproduction data. This leads to oj/^tv ~ 1 m b for E inc ~ 20 
GeV. However, the application of the VDM leads to a paradox |J - one obtains a to t{^'- 
N) « 0.7 • a to4 (J /ty-N), although, on the other hand, rv rs 2rj/*. This clearly indicates 
that the charmonium states produced in photoproduction are in a smaller - than average - 
configuration. Therefore, the VDM significantly underestimates a tot (J/^-N) |§. 

Indeed, the A-dependence of the J/\l/ production studied at SLAC at E inc ~ 20 GeV 
exhibits a significant absorption effect || leading to a a b s (J /^-N) = 3.5 ± 0.8 mb. It was 
demonstrated in [12| that, in the kinematic region at SLAC, the color coherence effects are 



still small on the internucleon scale for the formation of J/\J/'s and lead only to a small 
increase of the value of a a b s (J/^-N). So, in contrast to the findings at higher energies, at 
intermediate energies this process measures the genuine J/^-N interaction cross section at 
energies of ~ 15-20 GeV 0. 

To evaluate the nonperturbative QCD contribution we use an interpolation formula for 



the dependence of the cross section on the transverse size b of a quark-gluon configuration []. 
Three reference points are used to fix our parametrization of the cross sections (cf. Tab. |): 

• based on the observation that the <p-N cross section is nearly half of the p-N cross 
section we impose a(b ■ m p /m^) =10-12 mb, where bo is 0.6 fm. 

• 60 is approximately the transverse size of a tt, so we choose a(bo) = 25 mb, as known 
from 7T-N collisions. 

• for configurations with b > 1 fm we set cr(lfm) = 40 mb, the value reached when the 
two constituent quarks split to form two open charm mesons. 

Thus, the X-N cross sections is calculated via: 

a = J a(b) ■ \^(x,y, z)\ 2 dxdydz . (1) 

where y,z) is the charmonium wave function. In our calculations we use the wave 
functions from a non-relativistic charmonium model with a Cornell confining potential, 



V(r) = — - + see [I3| and refs. therein. 

The calculated cross sections of the charmonium states are shown in Tab. |. For the J/\I/ 
we found crj/^v = 3.6 mb. This is in good agreement with the data HJ, a — 3.5 ± 0.8 mb, 
for the J/^-photoproduction! 

Note the factor of ~ 2 between the < b 2 > values of the X10 versus the Xn, respectively: 
the transverse size of a fast x depends drastically on the polarization of the x states. Follow- 
ing the above discussion one therefore can expect significant fluctuations of the strength of 
the interaction due to this geometrically directed filtering, even in those kinematic regions 
where the color coherent phenomena seem to be at first sight only a small correction. The 
values of the x~N cross sections are given for different magnetic quantum numbers. 



1 Note that the deviation of a(b) from the estimate of PQCD occurs for b > 0.3 fm which is 
the distance close to the Chiral Phase Transition scale. We thank J. Bjorken for emphasizing this 
point. 



We calculate a(X-N) using the assumption that PQCD is applicable in this case. Here we 
will ignore the differences between bare quarks of the QCD Lagrangian and the constituent 
quarks. Then the expression for a(X-N) follows unambigously from the QCD factorization 
theorem [ ID| . The numerical calculation shown in Tab. | yields at CERN energies a(X- 
N) values which are significantly smaller than in 0. The <Jhard values Q in Tab. | are also 
calculated with eq. [I], but here we integrate only over the region 0.1 fm < b < 0.2 fm, because 
a s increases with b. Here b is the transverse distance between the heavy quarks transverse to 
the momentum of the heavy quarkonium. Thus, the region b > 0.2 fm cannot be calculated 
within PQCD. We choose the lower limit for the integral, b = 0.1 fm, because the quarkonium 
states will have a finite size already at the production point. For the calculation of the hard 
cross section the a (b) in eq. [I] is given by the function [[10 



a(b) = ^-b 2 a s (Q 2 )-xG(x,Q 2 ) . (2) 

The Bjorken x needed for the calculation of ahard(X — N) is calculated by x = 2 m N v > 
where Q 2 = ^, is the nucleon mass and v is the energy of the heavy Quarkonium state 
X in the rest frame of nucleon. The calculation in Tab. |I| was done for a state X produced in 
midrapidity at SPS energy and in the target fragmentation region at RHIC and LHC. One 
can see that the hard contribution to the cross section is just a correction at SPS energies, 
but at RHIC energies both contributions become comparable and at LHC it dominates (we 
neglect here that the DGLAP equation (Dokshitser-Gribov-Lipatov-Altarelli-Parisi) should 
be probably violated [|TT| ). 

We follow the analysis of (and refs. therein) to evaluate the fraction of J/\l/'s (in 
pp collisions) that come from the decays of the % an d No experimental information 
is available on the contribution from the x(1 + ~)j nor are an y contributions of the higher 



2 The calculated ahard reflects the effective cross section for the interaction of a QQ pair with 
nucleons of the residual nucleus, which will transform into the corresponding QQ state only after 
it has traversed the nucleus. 
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cc states (D-wave, radial excitations) known. Thus, the fraction of the directly produced 
J/^'s used in eq. |3] is an upper limit for the actual fraction of directly produced J/\l/'s. 
Furthermore, the ratio of directly produced J/^'s to the yield of J/\&'s due to \&' decays is 
J/ 1 !' : ^' = 0.92 : 0.08 in pp collisions 0. So, the suppression factor S of J/^'s produced 
in the nuclear medium is calculated as: 

S = 0.6 -(0.92- S J/ * + 0.08- 5*') + 0.4- . (3) 

Here are the respective suppression factors of the different pure charmonium states X 
in nuclear matter. The S x are for minimum bias pA collisions within the semiclassical 
approximation (cf. Hl5| ): 

% = A^N^X) =\j#B4, f{ B,.).^{-£ nf{ B,<W) ■ (4) 

Here p(B, z) is the local nuclear groundstate density (we used the standard parametrization 
from ||). 

The charmonium states are produced as small (r init « 0.2 fm) configurations predomi- 
nantly in gluon-gluon-fusion, then they evolve to their full size. Therefore, if the formation 
length of the charmonium states, If, becomes larger than the average internucleon distance 
(If > ttvat ~ 1-8 fm), one has to take into account the evolution of the cross sections with 



the distance from the production point [12 



The formation length of the J/\l/ is given by the energy denominator lj 



2p 



where p is the momentum of the J/^ in the rest frame of the target. With p = 30 GeV, 
the momentum of a J/\l> produced at midrapidity at SPS energies (E lab = 200 AGeV), this 
yields If m 3 fm. Due to the lack of better knowledge, we use the same If ~ 3 fm for the x- 
For the ^' we use If ~ 6 fm, i.e. we introduced an additional factor of 2, because it is not 
a small object, but has the size of a normal hadron, i.e. the pion. For = 800 AGeV we 
get another factor of two for the formation lengths due to the larger Lorentz factor. 

However, this has a large impact on the ^' to J/\l/-ratio depicted in Fig. [I], which shows 
the ratio 0.019 • SV/Sj/* calculated with (squares (200 GeV) and triangles (800 GeV)) and 
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without (crosses) expansion. The factor 0.019 is the measured value in pp collisions, because 
the experiments do not measure the calculated value Sw/Sj/y but J^^jy^ ■ B^ are the 
branching ratios for J/^, \&' — ► The calculations which take into account the expansion 
of small wave packages show better agreement with the data (circles) (taken from |2j) than 
the calculation without expansion time, i.e. with immediate J/^ formation, If = 0. We 
calculated this effect both at = 200 AGeV and 800 AGeV. The data have been measured 
at different energies {E lab = 200, 300, 400, 450, 800 GeV and ^ = 63 GeV). 

The charmonium survival probability S^ B in minimum bias AB collisions can be calcu- 
lated by |L5| S*b — &a ■ Sb- an d $b are defined by eq. (|. The neglect of the stopping 



of nucleons due to the energy loss seems safe: Drell-Yan processes do not suffer from this 
energy loss (see [16|,|l7j and refs. therein). 

The strong spin dependence of the x~N cross section is due to the angular dependence of 
the wave functions, which leads to different transverse quark separation for different states. 
The S-states (J/^/ and ^>') in contrast, do not exhibit any dependence on the angles. 

However, the P-states yield two vastly different cross sections (see Tab. |) for xio an d Xn> 
respectively. This leads to a higher absorption rate of the Xn as compared to the Xio- This 
new form of color filtering is predicted also for the corresponding states of other hadrons; 
e.g. for the bottomium states which are proposed as contrast signals to the J/\l/'s at RHIC 
and LHC! (For a detailed review of color coherence effects see [@.) Fig. ||a depicts the A 
dependence of the suppression of the Xn an d Xio i n pA an d in minimum bias AB collisions. 

The polarization of the x shown in Fig. |2|a may be difficult to observe. However, it is 
manifested in the production of all P-states (except in the J PC = ++ state, which in any 
case has a very small x ~^ J/^ + 7 branching ratio). Significant effects are expected for 
both 1 ++ and 2 ++ states, which give the major contribution to the J/\l/ production, and 
also for the S=0 state 1 H . This state was observed to date only in pp collisions ||19|| . The 
predicted effect is present for the 1 ++ state, since the probability of the m=0 state is 1/2 
for helicity +1 and is for helicity 0. 

There are other promising x states with the total angular momentum J=2. Fig. |2|b shows 



the different suppressions for these ^-states, i.e. with J=2 and J z = 0, 1(— 1) and 2(— 2). 
These states are mixtures of Xw-> Xn an d Xi-i'- x(Jj JJ = a i ' Xio + a 2 • Xn + a 3 ' Xi-i- The 
a's are given by the Clebsch-Gordan coefficients. 

The measured J /^/-suppression up to sulphur- uranium collisions seems to be described 
reasonably well with the present simple model, in view of the neglect of comoving 
hadrons ||[|. Note that the physically most reasonable scenario, namely the one including 
the expansion actually underestimates the suppression for SU by about 10 %. And also the 



apparently strong J/\l/ suppression found at CERN in going from SU to PbPb collisions J2G 
can not be explained in our calculations. This leaves space for more physics. Density fluctu- 
ations and the energy loss due to gluon radiation have also been neglected. Higher statistics 
for the pA-d&ta, are needed to clarify the importance of the expansion effects. For example 
the ratio of ^"s to J/\l/'s should be measured at one energy (e.g. y/s = 20 GeV, because 
here the Lorentz factor is not too large) and for different systems. This can prove that this 
ratio is not constant. One should also measure the rapidity distribution of this ratio to see 
the dependence on the Lorentz factor. 

We recommend to measure the A-dependent polarization of the \ a l so m pA reactions at 
•y/s = 20 GeV and in the midrapidity region, because the above described expansion effects 
and the interaction with secondaries will decrease this polarization. The polarization, once 
found, can be used to measure the interaction of these charmonium states with secondaries 
in AB collisions. The quantitative evaluation of the interaction with comovers, nucleon 
correlations and density fluctuations is postponed to the future. 

The other new effect predicted by the QCD factorization theorem is the fast increase 
of the cross section of spatially small configurations with energy, see Tab. |. Thus, QCD 
predicts a large nuclear absorption of charmonium in AB collisions at the RHIC and, in 
particular, at LHC energies, although the effective size of cc pairs will be ~ ^- ~ 0.1 fm 
due to the strong increase of the parton distribution functions at small x. The absorption 
depends also on the rapidity of the produced cc pair. 
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cc/66-state 


J/* 




AciU 


Acll 


T 




AolU 


Vli 1 

Aoii 


< b 2 > (fm 2 ) 


0.094 


0.385 


0.147 


0.293 


0.027 


0.15 


0.059 


0.12 


&nonperturbative (nib) 


3.62 


20.0 


6.82 


15.9 


0.31 


7.35 


1.48 


4.48 


a hard (mb) (SPS) 


0.024 


0.012 


0.021 


0.006 


0.101 


0.030 


0.103 


0.064 


Chard (mb) (RHIC) 


1.73 


0.68 


1.23 


0.30 


2.30 


0.64 


2.17 


1.27 


(mb) (LHC) 


20.8 


8.2 


14.7 


3.5 


28.4 


7.8 


26.3 


15.0 



TABLE I. The average square of the transverse distances of the charmonium states and the 
total quarkonium-nucleon cross sections a. For the x t wo values arise, due to the spin dependent 
wave functions (Im = 10,11). ahard are the perturbative QCD contribution at different energies 
(see text). 
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FIG. 1. The ratio 0.019 • 5V/Sj/<i> is shown in pA with the calculated cross sections for the J/^, 
Hf' and x (crosses) in comparison to the data (circles). The squares and the triangles shows the 
ratio calculated with the expansion of small wave packages (see text). 
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FIG. 2. The survival probabilities S for the different x-states are shown for pA and minimum 
bias AB collisions, a) gives the xio and the Xii(-i)j which differ by 0.2-0.3. The fraction of the 
Xio's of all x's increases to more than 50 % for P6P6-collisions. b) the x-states with J=2 and 
3 Z = 0, 1(— 1) and 2(— 2) are shown. 
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